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■ Abstract 

04 ■ A mixture of relativistic gases of non-disparate rest masses in a Schwarzschild metric is studied 

on the basis of a relativistic Boltzmann equation in the presence of gravitational fields. A BGK-type 

1 ^ ' , model equation of the collision operator of the Boltzmann equation is used in order to compute the 
^jH" non-equilibrium distribution functions by the Chapman-Enskog method. The main focus of this work 

I 1 is to obtain Fick's law without the thermal-diffusion cross-effect. Fick's law has four contributions, 

j^jfj. two of them are the usual terms proportional to the gradients of concentration and pressure. The other 

two are of the same nature as those which appears in Fourier's law in the presence of gravitational 
fields and are related with an acceleration and gravitational potential gradient, but unlike Fourier's 
law these two last terms are of non-relativistic order. Furthermore, it is shown that the coefficients 
of diffusion depend on the gravitational potential and they become larger than those in the absence 
of it. 
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1 Introduction 



The research on non-equilibrium properties of relativistic gases by using the Boltzmann equation in 
gravitational fields is a subject few explored in the literature. An important contribution was due to 
Bernstein [T], who obtained the constitutive equation for the non-equilibrium pressure of a relativistic 
gas and the corresponding transport coefficient of bulk viscosity in a Friedmann- Robertson- Walker metric. 

Recently a relativistic gas in a gravitational field was analyzed in order to determine the influence 
of the gravitational potential gradient on Fourier's law [2] [3J and of the influence of the gravitational 
potential on the transport coefficients [3j. 

According to [3J the heat flux in Fourier's law has three contributions, the usual temperature gradient 
and two relativistic terms. One of them - proposed by Eckart from a thermodynamic theory [1] - is 
connected with the inertia of energy and represents an isothermal heat flux when matter is accelerated. 
The other - suggested by Tolman [SJ [B] - requires that in absence of an acceleration a state of equilibrium 
of a relativistic gas in a gravitational field is achieved only if the temperature gradient is counterbalanced 
by a gravitational potential gradient. 

An open question that follows refers to the modification of Fick's law in the presence of gravitational 
fields and the answer is the aim of this present work. 

In this work we are interested only in analyzing Fick's law in the presence of gravitational fields 
without the determination of the thermal-diffusion cross-effect and of the heat flux with the corresponding 
diffusion-thermal cross-effect, which will be subject of a forthcoming paper. As in the work [3J we use 
the Schwarzschild metric and analyze a relativistic gas mixture of constituents which have non-disparate 
rest masses. 
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Here we show that Fick's law in the presence of gravitational fields has four contributions. The usual 
contributions due to the concentration gradient and pressure gradient, apart from the two that appear 
in Fourier's law and which are proportional to the acceleration and gravitational potential gradient. 
However, unlike Fourier's law these two last contributions are not of rclativistic order. We have also 
obtained that the diffusion coefficients depend on the gravitational potential, becoming larger than the 
ones in the absence of it. 

The work is structured as follows. In section 2 we introduce the Schwarzschild metric, the system of 
Boltzmann equations in the presence of a gravitational fields and the two first moments of the distribution 
functions with their corresponding balance equations. A BGK-type model of the Boltzmann equation 
is introduced in section 3, which depends on a reference distribution function, determined from the 
assumption that the balance equations of the full Boltzmann equation and of the BGK-type model lead 
to the same production terms. The non-equilibrium distribution functions are calculated in section 4 
by using the Chapman-Enskog method. In section 5 the constitutive equations for the diffusion fluxes - 
which correspond to Fick's law - are determined from the non-equilibrium distribution functions and the 
diffusion coefficients are obtained. In the last section the main conclusions of the work are stated. We 
close the work with two appendices. In the first one it is shown how to calculate the production terms 
of the partial balance equations of the energy-momentum tensors, while in the second the components of 
the Christofcll symbols in a Schwarzschild metric are given. 



2 Basic Equations 

We consider a relativistic gas mixture of r constituents in a Ricmannian space with metric tensor g^ 
where the particles of constituent a = 1, . . . r have rest mass m a and are characterized by the space-time 
coordinates = (ct, x) and momenta (p£) = (p°,p ). The length of the momentum four-vector is 
9fj,vPaPa — m a c2 i which implies that 



PaO 90iPa / ■?•?,/ \ 4 j /1 \ 

Pa = 1 PaO = V 900 m c + \90i90j ~ 9009ij)PaPa- (1) 

5oo 

As in the work [3] we shall adopt the isotropic Schwarzschild metric 

GM 
2c 2 r 



ds 2 = g (r) (dx Q f g 1 (r)6 ij dx i d^, g (r) = V fg 2 , 9l (r) = jl+ , (2) 



where G denotes the gravitational constant and M the mass of the spherical source. 
In terms of the isotropic Schwarzschild metric @, eqs. ([T]) reduces to 



Pa = — PaO = \/3oV TO a c2 +ffl|p a | 2 , V^fl = J 9o9\- (3) 

go v 

where g = det{g^ v ). 

The components of the four- velocity in the isotropic Schwarzschild metric are 

Jjn = ( — — C V | (4) 

V V9o ~ v 2 /c 2 ' a/90 ~ v 2 /c 2 J ' 

which in a comoving frame (v = 0) reduces to = (c/^/go, 0). 

A state of the relativistic mixture of r constituents in the phase space spanned by the space-time 
and momentum coordinates is described by the set of one-particle distribution functions /(x, p a ,t) = 
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f a , (a = 1, 2, . . . , r) such that /(x, p Q , t)d 3 xd 3 p a gives at time t, the number of particles of constituent a 
in the volume element d 3 x about x and with momenta in the range d 3 p a about p a . 

In the presence of gravitational fields the one-particle distribution function of constituent a satisfies 
the Boltzmann equation (see e.g. [7]) 

Pa^A ~ VaPajk = E / (/»/& 1 ah)F ha O ah di}^^ . (5) 



Here r* are the Christoffel symbols, Fb a = (paPb/iY — '"- a " l b l 
a a b and dQ denote the invariant differential elastic cross-section and the element of solid angle that char- 
acterizes a binary collision between the particles of constituent a with those of constituent b, respectively. 
The binary collision is characterized by the momentum four- vectors of the particles of two constituents 
p" and p" before collision and p'^ and p' b a after collision, which obey the energy-momentum conservation 
law Pa+Pb= p'a + v'b ■ Furthermore, the following abbreviations were introduced in ([5]): 

f' a = f{fr V ' a ,t), f b = /(X,PU), fa = f(*,Pa,t), / 6 = /(x,p 6 ,f). (6) 

The two first moments of the distribution function are the partial particle four-flow N£ and the partial 
energy-momentum tensor Tg v of constituent a, which are defined in terms of the one-particle distribution 
function by 

K = cf rtf ay /=gipL, rr = c [ pX/.>/=s— • (7) 

J PaO J PaO 

By summing the partial quantities we get the particle four-flow of the mixture iV M = J2a=i ^a an d the 
energy-momentum tensor of the mixture T Ml/ = Xm=i T£ v . 

The balance equations for the partial particle four-flow and for the partial energy-momentum ten- 
sor of constituent a are obtained by multiplying the Boltzmann equation ([5]) by Cy/—gd 3 p a /p a o and 
cpay/—gd 3 p a /pao, and integrating the resulting equations. Hence, it follows 



a \fi 

6=1 



0, Tr. v = E C / (P'a - PZ)fafbFba*a b dn^g-^V=9—, (8) 

J PbQ Pan 



where the semicolon denotes the covariant derivative. 

The balance equations for the particle four-flow and energy-momentum tensor of the mixture are 
obtained through the sum of (jHJ over all constituents, namely, N 1 *-^ = and '. v = 0. 

In this work we are interested in analyzing the diffusion of a relativistic gas mixture in the presence 
of gravitational fields. For that end we shall derive only the constitutive equation for the diffusion flux, 
which corresponds to Fick's law. The diffusion flux of the constituent a in the mixture is introduced 
through the decomposition of the the particle four-flow as 

N% = n a C/'' + J£, such that J£[/ M = 0, (9) 

where n a denotes the particle number density and with U^U^ = c 2 is the four- velocity. In terms of 
the projector 

A ^ = gilv _ jLumr with a^[/„ = o, (10) 

the diffusion flux can be represented - thanks to (JTJ) i - as 

% = [tifa — . (11) 

J PaO 
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The sum of © over all constituents by taking into account the decomposition of the particle four-flow 
of the mixture TV' 1 = n£/ M - with n denoting the particle number density of the mixture - implies that 



n = E n - E J « = °- ( 12 ) 



a=l 



Hence, due to the constraint (IT21) 9 there exist only (r — 1) partial diffusion fluxes that are linearly 
independent for a mixture of r constituents. 



3 BGK-type kinetic model 

In kinetic theory of gases it is usual to replace the collision operator of the Boltzmann equation by models 
which arc more easy to handle. These model equations are referred generally as BGK-type kinetic models 
and the first one was proposed independently in 1954 by Bhatnagar, Gross and Krook [8] and by Welandcr 
[3]. The extension of BGK-type models to relativistic gases was due to Marie [TUJ[TT] and Anderson and 
Witting [H]. 

The main objective of this work is to discuss the effect of gravitational fields in the diffusion of 
relativistic gas mixtures and to derive Fick's law in the presence of gravitational fields without the cross 
effect of thermal-diffusion. For that end we shall use a truncated Grad's distribution function that 
depends only on the diffusion flux, which will be the only non-equilibrium quantity in this theory. The 
consideration of the full Grad's distribution function is more involved and this will be the subject of a 
forthcoming paper, where the heat flux and cross-efects shall be analyzed. 

The model equation we shall use is based on the McCormac model for mixtures of non-relativistic gas 
mixtures [13] extended to a relativistic truncated Grad's distribution function. The McCormac model for 
non-relativistic gas mixtures of non-reacting and reacting gases by using a truncated Grad's distribution 
function were analyzed in the works |14j and |15| . respectively. 

The truncated Grad's distribution function which takes into account the diffusion flux as the only 
non-equilibrium quantity is given by (see [TB]) 



h = fi 0) (l-^Pa,) 
V Pa / 



(13) 



where p a = n a kT is the pressure of constituent a and fd the Maxwell- Jiittner distribution function 



f (0) = n JL exv (- l M.\ = ^ cxpf c V m a c2 +9M 2 ) (u) 

Ja A-KkTmlcK 2 {Ca) P \ kT J A-KkTmlcK 2 {C a ) P \ kT J' y ' 

Here the second equality is written in a comoving frame, where = (c/^go, 0). Above T denotes 
the temperature of the mixture and it was supposed that all constituents are at the same temperature. 
Furthermore, k is the Boltzmann constant and Ca = m a c 2 /kT represents the ratio of the rest energy 
m a c 2 of a relativistic a-particle and the thermal energy of the mixture kT. The limiting case of low 
temperatures £ a 3> 1 corresponds to a non-relativistic regime and the limiting case of high temperatures 
(„ < 1 to an ultra-relativistic regime. Furthermore, K n (Ca) denotes the modified Bessel function of 
second kind which is given by 

K n (C a )=(^T^L f e-W*/ 2 -ir 1/2 ^ (15) 



2 / r(n + l/2) 
The BGK-type kinetic model we shall use reads 



Pa^~ r >>«H = " 2 ^ ^ ~ ^ (16) 
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where m a b = m a mb/(m a + rub) is the reduced mass, v a b is a collision frequency whose inverse is of order 
of the mean free time and f^ b is a reference distribution function. In order to be compatible with the 
truncated Grad's distribution function the reference distribution function must be written as 



(17) 



where B a b and C^ b are coefficients to be determined. 

In order to determine the coefficients B a b and C^ b we multiply the BGK-type kinetic model ([TBI) by 
c^Z—gcfipa/pao and cp v a yj— gd 3 p a /p a o and integrate the resulting equations. Hence, it follows 



K-^ = -2 E ^ ">* C fifa- fab)\ 
6=1 J 

6=1 J 



d 3 Pa 
PaO ' 



-d 3 p a 
1 • 

PaO 



(18) 
(19) 



Comparison of the above equations with © leads to 



- 2 V m ab v ab c I \f a - fabjV^g—^- = 0, 
r-i. J Pao 

6=1 



2 m ab Vab C f P%(fa ~ /3)a/=0— = £ C / (f# - f£) J 'afbFbaO ^dSlyf^ 
6=1 J Pa0 6=1 J 



d 3 p b , d 3 p a 

9- 



PbO 



PaO 



(20) 
(21) 



and from the above equations we can determine the coefficients B a b and C" h . Let us denote the integrals 
in igOJl and dUJ as 



3 T M 
x a6 



C (fa- f*b)V~9 



d 3 p a 
PaO 



2 Kb = clpZ(fa-fDV—g 



d 3 p a 
PaO 



C I (PT -p^fafbFbaOabdnV—g^V—g-^ 



PbO 



PaO 



The integrations of lead to Q 

IT - Pa ( B ab ^l(Ca) TO a ^ \ 

ab - m a \kT K 2 (Ca) + 

K 3 (Ca) m a 



X a6 Pa 



B a 



kT~ ' K 2 ((a) kT LabU " U p a 



■17" 



(22) 
(23) 

(24) 
(25) 



thanks to ([T3|) and (|17|) . The integral (|23[) is calculated in the Appendix A for mixtures of non-disparate 
masses (mj w m a ) and constant differential cross section cr a b. Its expression is given by 



/ a 6 

x n6 — ± a6 I 



J a _ ^ 
Pa P& 



where the following abbreviation was introduced 
a 167rp a p b cr ab 



3cA- 2 (Ca)i^2(C6) 



-J=^3 (2VU6) + ( 2 + 7V ) ^2 (2 VCaO 
VCaC& V ' V CaC&/ V 



(26) 



(27) 



1 For more details on the calculations of these integrals one is referred to the work [3]. 
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From 



(|2"B)) we conclude that C^ b U^ = 0, B ab = and 



°a6 



T 



a b 



I 1 ' 
J b 



2m ab v ab p a \p a p b 
Hence, the reference distribution function (|17[) becomes 



f R = f (0) 

J ab J a 



_2m ab v ab p a \p a p b 



I'' 
J b 



J'' 
Pa' 



J'' 
-'a 

Pa 



Par 



(28) 



(29) 



4 Non-equilibrium distribution function 

For the determination of the non-equilibrium distribution function we shall rely on the Chapman-Enskog 
method. First we write the one-particle distribution function of constituent a as f a = /„ (1 + fa) where 
the deviation of the Maxwell- Jiittncr distribution function is considered a small quantity (\f a \ < !)• 
According to the Chapman-Enskog method the equilibrium distribution function fi°^ is inserted on the 
left-hand side of the BGK-type kinetic model (|16p and the representation f a = /i ' (1 + ip a ) on its 
right-hand side. Hence, by performing the derivatives we get that the deviation of the Maxwell- Jiittner 
distribution function becomes 



1 



fa 



Pa dn a 

2 121=1 m <>b "ab I n a dx u 

c 2 dgip a p a p k a x l 

OjiOhl — 

2kT dr U T p T a J r 



1 - 



K 3 (Ca)Ca , PlUr 
K2((a) 



__ 

T 

J_ c r i pIpIp" 

k T 9l0ij ™ U T pl 



kT 



2JmkT 



8T plpl dU T 
dx v kT dx u 

PaPa— S 



dr 



'ij 



6=1 



Pa 



IP IP 

__ _ zJl I 

Pa Pfa 



K 

Vaii- 

Pa 



(30) 



Note that the deviation of the Maxwell- Jiittner distribution function (p a is a function of the diffusion flux 
J^, gradients of n a , U^, T and derivatives of the components of the metric tensor go, g\. 



5 Fick's law 

Insertion of the representation f a = fa°^ (1 + fa) together with (TT4")) and (ptUl) into the definition of the 
diffusion flux and integrating the resulting equation it follows a relationship which can be used to 
determine Fick's law, namely, 



6=1 



/ If" |M 

J2 T * b ( — - — 

\Pa P6 



dp a _ n a h a 

dx v c 2 



A^U T U U . T 



1 <9$ 
_ A M — 

1 - $ 2 /4c 4 dxP 



(31) 



Here h a = m a c 2 K^(C, a ) / .K^CCa) denotes the enthalpy per particle of constituent a and <E> = —GM/r the 
gravitational potential. 

If we sum (f3"T|) over all constituents we get the momentum density balance equation for the mixture 
(see [3] eq.(22)) 



A 



nh 

72" 



ts v u T u v . r - 



1 <9$ 
_ A w — 

l-$ 2 /4c 4 9a;J 



= 0, 



dx u 

where p = 12 a= i P a i s the pressure and h = 12 a —i n ah Q /n the enthalpy per particle of the mixture. 



(32) 



(i 



Since we have r — 1 diffusion fluxes which are linearly independents we rewrite (|3ip as a system of 
r — 1 equations for the determination of the diffusion fluxes: 



6=1 



0a" + lXa j ftc" 



n a h a - nh 



c 2 p 



A' ll/ [/ T £/„. T 



1 



(33) 



Here we have used the constraint (|32[) and introduced the concentration x a = p a /p 
a and the coefficient 



D 



ah 



In a comoving frame where 
A 00 = 



A 11 



Sab 
Pa 



Pr 



Pfe 



T 



.g*3 



we obtain from (1331 Fick's law: 



j> 



r-1 

£ 

6=1 



nZ> o6 <J* 



+ (x 6 - 1) 



<9 hip 



.9i 



nth 



nh 



c 2 p 



\/ 2c 2 



U T U VT 



-S ij 

I u 1 



n a /n of constituent 
(34) 

(35) 



<9$ 



1 - $ 2 /4c 4 dxi 



(36) 



Here we identify V a b = D a b/(ngi) as the diffusion coefficient. 

From Fick's law (|36[) we may infer that the diffusion flux - without the consideration of thermal- 
diffusion effects - has four different contributions: 

1. a flux due to a concentration gradient that tends to reduce the non- homogeneity of the mixture; 

2. a flux associated with the pressure gradient where heavy particles tend to diffuse to places with 
large pressures as in the case of centrifuges; 

3. a flux related with an acceleration term, since the acceleration acts on particles with different 
masses; 

4. a flux caused by the gravitational potential gradient. By taking into account that the acceleration 
term is absent and that the pressure is uniform, we may infer from (|36|) that the diffusion flux 
vanishes if the concentration gradient is counterbalanced with the gravitational potential gradient. 

Let us compare Fick's law with Fourier's law for a single relativistic fluid. According to cq. (33) of 
[3] we have 



J dT T 
\ dx J c 2 



U T U VT - 



l_$2/ 4c 4 9a .j 



(37) 



where q l denotes the heat flux and A the thermal conductivity coefficient. By comparing (|3"6"|) with (|3T|) 
we infer that the combination of an accelerated term and a gravitational potential gradient has the same 
structure in both equations. However, in the case of the heat flux this combination is of relativistic order 
which is not case of the diffusion flux, since in the non- relativistic limiting case (£ a ^> 1) the expression 
(n h — nh)/c 2 p tends to (x m — X) a =i x am a )/kT. As was pointed out in [3] the acceleration and the 
gravitational potential gradient terms, which appear in the heat flux (|37|). were first proposed by Eckart 
[4] and Tolman (5J [6] , respectively. 
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In a non-relativistic theory (see [T7] pp. 344 and 345) the diffusion flux is proportional to the so-called 
generalized diffusion forces whose expression, in our notation, is given by 



m„n„ 



Eb TO b n « 



Vlnp - 



J2 b m b n b F t 



(38) 



Here, F a denotes a force per unit mass that acts on the particles of constituent a. Both equations 
and (p?5j) depend on the concentration and pressure gradients, but the third term of ([55)) is of different 
nature. Indeed, it is connected with an acceleration when the forces acting on the particles of different 
constituents are unequal. However, this term vanishes when one considers that only gravitational forces 
are acting on the particles. 

One important point to call attention is that the acceleration term is always connected with the 
gravitational potential gradient. If we eliminate from (|3^|) the acceleration term by the use of the 
momentum density balance equation for the mixture (|32[) both contributions disappears and we get 



r-1 

E 

b=l 



nbhb^ Ship 



nh 



(39) 



Note that the terms in brackets reduce to the first two terms of (|38|) in the non-relativistic limiting case. 
Although this is a common practice in the non-relativistic and special relativistic kinetic theories, here 
this operation cannot be done, since in this case the dependence on the gravitational potential gradient 
disappears from the constitutive equation of the diffusion flux and in our opinion this is not physically 
correct. This fact was also shown in [3J, that the elimination of the acceleration term from Fourier's law 
leads to a heat flux which is independent on the gravitational potential gradient. 

As in the case of a single relativistic gas in the presence of gravitational fields [3] the diffusion coefficient 
depends on the gravitational potential through the component of the metric tensor gi(r) = (1 — $/2c 2 ) 4 . 
Let us analyze the diffusion coefficient of a binary mixture which follows from (|34[) . In this case we have 
that 



P1P2 



3cX 2 (Ci)if 2 (C2) 



nTi 2 167rncri2 (1 - |$|/2c 2 ) 4 
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VC1C2 



C1C2 



(40) 



From the above equation we infer that the diffusion coefficient becomes larger in the presence of gravita- 
tional fields. As was pointed out in [3] this increase is very small for stellar objects like Earth and Sun, 
but it becomes important for more compact stellar objects, like a neutron star. 

The diffusion coefficient (14TJ1) in the non-relativistic limiting case (1,(2 3> 1 that corresponds to a 
binary mixture of constituent with non-disparate masses (mi w 7712) at low temperatures and under the 
influence of weak gravitational fields reads 



V 



11 



32ncri2 



kT 



1 - 



11 

I6C1 



1 



2|$| 



(41) 



Without the relativistic and gravitational field corrections given by the underlined terms of (|4ip we get 
the self-diffusion coefficient of a mixture of hard sphere particles, since in this case the differential cross 
section is given in terms of the relative diameter of the particles at collision by a 12 = {di + d2) 2 /16 (see 

e.g. nanFj). 

For high temperatures £ Q <C 1 and the diffusion coefficient (j40]l in the ultra-relativistic limiting case 
under the influence of weak gravitational fields becomes 



107rneri2 



1 



10 



1 



2|$l 



(42) 
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6 Conclusions 



The main objective of this work was to show that Fick's law in the presence of gravitational fields has 
the same structure as Fourier's law for a single constituent, i.e., apart from the well known forces without 
the presence of the gravitational field, both laws depend on the acceleration and gravitational potential 
gradient. From the analysis of Fick's law we have concluded: (i) the driving forces in the absence of 
gravitational fields are the usual terms connected with the concentration gradient and pressure gradient; 
(ii) the terms related with the acceleration and gravitational potential gradient unlike Fourier's law are 
of non-relativistic order; (iii) the transport coefficients of diffusion depend on the gravitational potential 
and they become larger than the ones in the absence of it. 

A Appendix: Evaluation of 3 Z^ 6 

In order to perform the integration of 



-d 3 p b ^r-^ d3 P a 

PbO PaO 



3 l» ab = c I (pT - p£)/ / 6 F fa CT a6 dlV^— %/=S— , (43) 



we introduce the total momentum four-vector PP and the relative momentum four- vector defined by 
(see e.g. [19]) 

P" = pS+p£, p'^pT+pT, Q»=rt-tf, Q"*=pT-pT- ( 44 ) 

The momentum four-vector conservation law when applied to the above equations lead to 

pf = P'^ P^Q^iml-mDc 2 , Q 2 = P 2 -2(ml + m 2 b )c 2 . (45) 

Here P 2 = P M P M and Q 2 = —Q^Q^ are the magnitudes of the total and relative momentum four-vectors, 
respectively. Furthermore, the inverse transformations of (|4"4"|) read 

Pa = — + —, Pb=^ 2"' Ptt 2 IT' Pb 2 2~' ( ' 

and the Jacobian of the transformation from (p%,p b ) to (P M , Q M ) is 1/8, so that d 3 p a d 3 p b = d 3 P d 3 Q/8. 

By introducing a space-like unit vector k M which is orthogonal to P M - i.e., k''P M = - the relative 
momentum four-vector can be written as 

QM = ( W a-^)c 2 pA1 + k^ ^ pA _ 2p2(m 2 ~ m 2 )c 2 + (m 2 _ m 2)2 c 4. (47) 

As in the work [16] we shall restrict ourselves to the case where the rest masses of the particles of the 
constituents are not too disparate so that mt « m a , (V6 = 1, . . . , r). In this case the relative momentum 
four- vector and its modulus can be approximated by 

Q^ = Q k", Q 2 =P 2 - 4m a m b , (48) 

Furthermore, in this approximation P^Q^ = 0, and the the invariant flux reduces to 

F ba = ^(pZ Pb ,) 2 - mlmy = (49) 

In terms of the total and relative momentum four- vectors the linearized product of the distribution 
functions can be written as 



, , n a n b e 

Jajl 



16^ 2 fc 2 T 2 C 2 m>2if 2 (Ca)^2(Cb) 



2 Pa 2 Pb 



(50) 
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thanks to {13} and 

For the integration of (|43[) with respect to the relative momentum four-vector, the center-of-mass 
system is chosen where the spatial components of the total momentum four- vector vanish. By introducing 
the representations (P^) = (P°,0) and (Q a ) = (0, Q) we can obtain: 



ha 



d 3 p a d 3 p b _ 1 d 3 Pd 3 Q 
Pao Pbo 4:^90 Po 



(51) 



Now we write the element of solid angle as d£l = sin 0^0^$, where and $ are polar angles of 
the spatial components of Q' a with respect to Q a and such that O represents the scattering angle. The 
differential cross section cr ab depend on the modulus of the relative momentum four-vector Q and on 
the scattering angle 0, i.e, a ab = cr ao (Q,&). Here we are interested in analyzing a mixture where the 
differential cross section is constant, which in the non-relativistic case corresponds to a mixture of hard 
spheres particles. Hence without loss of generality, we may assume that the spatial component of Q a is 
in the direction of the axis a; 3 , so that we can write Q a and Q' a as: 



(Qn = q 



o 



W 



(Q'n = q 



( 

sin 6 cos 4> 
sin sin $ 
V cos 9 



(52) 



(53) 



(54) 



where dfl* = sin 8d8d4> represents an element of solid angle. Hence, we can obtain the following expres- 
sions for the integrals 

r 2,T 



From the above considerations we can obtain the following result 

J (pjf -p£)dn =\J (Q ,fJ - - Q**)dQ = -2nQ». 
We can also perform the integrations in the spherical angles of Q^, denoted by 9 and 0, i.e. 

/ . 

(Q") = Q S ™i C ™t | , d 3 Q = Q 2 sm0d6d4>dQ = Q 2 dn*dQ, 




dn* 



sin 9d6d<j) = Air, 



Q^dVL* = 0, 



Q»Q v dW = — Q 2 



By using the above results the integral (|4"3"]) can be written as 



3 T M 



n a n b a a b 



^av ^bv 

48k 2 T 2 cm 2 a m 2 K 2 (( a )K 2 (C b ) V~p7 ~ ~Pb~ 



p»u„ (P^pv 

kT I 

V p 2 



.All' 



Q°dQJgf go 



-d A P 



(55) 



(56) 



For the integration of (|56|) we consider a Lorentz rest frame where (U^) = (c/y/go, 0), a spherical 
coordinate system and write d 3 P = |P| 2 d|P| sm$ dip d$. Furthermore, we introduce a new variable of 
integration 



KaCb 



Q 2 



m a mbC 



so that 



kT 

QdQ= ( — 



: dx 



(57) 



where the range of this new variable is given by 2-^/CaCf) < x < oo and write the time and spatial 
coordinates of the total momentum four-vector as 



Po 



kT 

c 



9o xy, 



kT x\J y 2 — 1 



9i 



(58) 
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In this case the element of integration becomes 



d 3 P ( kT\ 2 

9i9o-pT={—J x 2 y v / y 2 -l sin d dip d-d dy, 



(59) 



and the range of integration of the variable y reads 1 < y < oo. 

From the integration over the solid angle and over the variable y we can obtain the following results 



^ty d 3 P 



9i9o e kT — = 4tt — xK^x) 



Po 



IT 



c 



'9i9a I 

so that ([56]) reduces to 
An a n b a ab 



r»u„p»p»d 3 P fkT\ 



u»u v 

xK z {x)—z--K i (x)gr 



3tM 



3cK 2 (( a )K 2 (( b ) \p a p b JJ 2 



J a J b 



'CaCb 



[x 2 K 3 (x)-3xK 2 (x)] 



16C Q 2 C 6 2 2C a C b 



(60) 
(61) 

dx. (62) 



To perform the last integration in the variable x in (|62| we need the following integrals of Bessel 
functions 



poo p oo 

/ x 5 K 2 (x)dx = x 4 [xK 5 (x)-QK 4 ( X )}, / x 3 K 2 (x)dx = x 3 K 3 ( X ), 

I poo 

xK 2 {x) dx = X Ki{x) + 2K (x), / x e K 3 (x) dx = x 5 [xK e (x) - 8K 5 ( X )] , 



x 4 K 3 (x)dx = x A K i ( X ), 



x 2 K 3 (x) dx = x 2 K 2 ( X ) + 4x^i (x) + 8K (x)- 



Hence (l6"2l becomes 



37-/^ x [ a b 



ZJL _ _ 
Pa Pb 



where we have introduced the following abbreviation 



I6n a n b a ab 
3cK 2 (( a )K 2 (( b ) 



(63) 
(64) 
(65) 

(66) 
(67) 



B Appendix: Isotropic Schwarzschild metric 



For the isotropic Schwarzschild metric (see e.g. [2"0] ) 
ds 2 = 9o(r) (dx ) 2 - gi(r)S ij dx' l dx J , g (r) = 



1c l r 1 



GM ' 



. , / GMy 



the Cristoffel symbols read 



rE!o = 0, 
rf 

r 



r u - = 



r"=o (i^j jtk), rj, = o, 



1 dgi(»-) g ^_ 

" ~ 2g x (r) dr jk r ' 

_J dg (r) x i 

2gi(r) dr r ' 



r 



0; 



0j 

1 dg (r) x 3 
-Oi-, — 



00 



2g (r) dr 
1 dgi (r) x 3 
2gi(r) dr r 



(69) 
(70) 

(71) 
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where the underlined indices are not summed and 



dg (r) _ 2GM (1 - dgi(f)_ _ 2GM / GM\ 3 

dr " c 2 r 2 (i i GM\ 3 ' dr c 2 r 2 \ 2c 2 r J 

\ 2c' 2 r J 
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